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ABSTRACT 

The consortium for Galactic studies with the Arecibo L-band Feed Array (ALFA) is conducting a 
neutral hydrogen (HI) survey of the whole Arecibo sky (declination range from —1° to 38°), with high 
angular (3.5') and velocity resolution (0.2 km s~ x ). The precursor observations with ALFA of a region 
in the Galactic anti-center reveal numerous isolated, small (a few pc in size), and cold (X^ < 400 K) HI 
clouds at low negative velocities, distinctly separated from the HI disk emission ('low- velocity clouds', 
LVCs). These clouds are most likely located in the transition region between the Galactic disk and halo 
(at scale heights of 60-900 pc), yet they have properties of typical cold neutral clouds. LVCs are colder 
and, most likely, smaller and less massive than Lockman's clouds in the disk/halo interface region of the 
inner Galaxy. Our observations demonstrate that the cloudy structure of the interface region is most 
likely a general phenomenon, not restricted to the inner Galaxy. LVCs have sizes and radial velocities 
in agreement with the expectations for clouds formed in low-temperature fountain flows, although we 
measure a factor of ten higher HI column densities. Alternatively, LVCs could represent the final stages 
of the infalling intergalactic material in the on-going construction of the Galaxy. 

In the same dataset at higher negative velocities, we have discovered a 'companion' HI cloud located 
50' southwest of CHVC186+19-114. CHVC186+19-114 is a typical compact high velocity cloud (HVC) 
with a well-defined core/envelope structure. The companion cloud has a diameter of only 7' x 9', and 
is one of the smallest HVCs known, most likely stripped from the main cloud through the interactions 
with the halo medium. 
Subject headings: ISM: clouds 
medium 
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1. INTRODUCTION 



The Galactic halo, a hot, ubiquitous multi-phase gas, 
with a low density, and temperature around 10 6 K, that 
surrounds our Galaxy is crucial to Galactic activity. As 
well as acting as a thermal insulator, the halo constrains 
models of the interstellar medium (ISM) in the Galactic 
disk through its impact on the dynamics of the hot in- 
tcrcloud component. From the theoretical point of view, 
the Galactic fountain model (Shapiro & Field 1976) was 
the first to address interactions between the disk and the 
halo. In this model, the hot gas flows from the disk into 
the halo, starts to cool, and falls back onto the disk. One 
of the later models, so called chimney model of the ISM 
(Norman & Ikeuchi 1989), suggested that the disk and 
the halo are connected through chimneys - large conduits 
resulting from superbubbles bursting out of the disk and 
forming collimated structures that can reach heights of 
~ 1 kpc. From the observational point of view, numer- 
ous ISM structures (shells, superbubbles, worms, chim- 
neys, plumes) were found in the vertical gas distribution 
of our own Galaxy and several nearby galaxies, and have 
been postulated to play an important role in the transfer 



of mass and energy between galactic disks and their ha- 
los. However, the details of how these structures replenish 
the hot halo gas are still not fully understood (Dove et al. 
2000). 

In understanding the ways in which the Galactic disk 
and halo are related the transition region between the 
Galactic thin disk and the halo, also known as the Galac- 
tic thick disk, is of special interest. This region extends to 
heights up to 1—1.5 kpc above the plane. The disk/halo in- 
terface region was discovered several decades ago (Shane 
1967; Lindblad 1967) and was considered to represent a 
smooth envelope of neutral hydrogen (HI) surrounding the 
Galactic spiral structure. Further studies (Albert 1983; 
Danly 1989; Lockman 1986) found occasional isolated halo 
clouds at lower heights and suggested that this component 
could be pervasive, yet patchy, and most likely dynami- 
cally connected with the disk. 

Galactic fountain and chimney models, as well as recent 
numerical simulations (de Avillez 2000), predicted the ex- 
istence of small-scale, cloudy structure in the interface re- 
gion, with clouds ranging in size from a few parsecs to 
tens or hundreds of parsecs. Recent observations with the 
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Green Bank telescope (Lockman 2002) have actually re- 
vealed dramatic small-scale structure of the interface re- 
gion in the inner Galaxy, composed of numerous discrete 
HI clouds. While this confirms the expectations from the- 
oretical models and simulations, the full extent and the 
origin of this cloudy component is still not understood. 

Another set of attractive questions concerning the 
disk/halo interface reagion comes from the idea that the 
formation of galaxies is an ongoing process (Oort 1966) 
and that the accretion of intergalactic material (IGM) 
serves as the main source of new star formation fuel 
(Mailer & Bullock 2004; Murali ct al. 2002). However, 
the details of the transformation from the smooth and hot 
IGM on kpc scales, down to clumpy and cold gas on star- 
forming scales, still need to be worked out. For example, 
what is the internal structure of IGM clouds at different 
stages of their infall? As searches for accretion remnants 
around other galaxies are difficult due to sensitivity and 
area limitations (e.g., Fraternali et al. 2002), our own 
Galaxy is still one of the best places to start seeking an- 
swers to these questions. 

HI clouds in the Galactic halo at heights > 1 kpc have 
been studied observationally for several decades and are 
often classified into several distinct groups based on how 
much their velocities deviate from Galactic rotation (in- 
termediate and high velocity clouds, HVCs), and their 
relative size (compact and extended clouds). Velocities 
reach approximate +/ — 450 km s _1 in the LSR reference 
frame and their sizes vary from over 100 deg 2 to com- 
pact HVCs with sizes of 4-30 arcmin (Putman et al. 2002; 
Brims & Westmeier 2004). In addition, almost all HVCs, 
especially the more compact ones, clearly show head/tail 
or core/envelope structures, with heads or cores having a 
narrow linewidth while tails and envelopes having a signif- 
icantly broader velocity structure, possibly showing signs 
of interaction with the ambient medium (e.g. Brians & 
Mebold 2004). It is still not clear whether this range of ob- 
servational properties, especially the low-end of the cloud 
size and velocity distributions, is caused by observational 
limitations. What determines the size and morphology of 
the observed halo clouds? Are these remnants of very dif- 
ferent physical processes, or do different 'types' of clouds 
represent different stages of the same phenomenon — ac- 
cretion onto the Galaxy, for example? 

The questions above are just a subset of issues that still 
await to be addressed. The main purpose of this paper 
is to summarize the ongoing observational efforts to un- 
derstand the properties of gas clouds from the halo to the 
plane by studying the disk/halo interface of our Galaxy. 
This work is inspired and enabled by the new multi-beam 
system, very recently installed at the Arecibo telescope 8 , 
which allows sensitive, large-scale HI surveys with an an- 
gular resolution of 3.5'. Previous large-scale surveys of 
the link between the disk and the halo have been limited 
by resolution (Lcidcn-Dwingcloo Survey, 36'; Hartmann & 
Burton (1997)), or an inability to trace features from the 
halo into Galactic emission (Putman ct al. 2002). Recent 
high-resolution Galactic Plane surveys (Taylor et al. 2003; 
McClurc-Griffiths et al. 2005) have been limited to trac- 

8 The Arecibo Observatory is part of the National Astronomy and 
agreement with the National Science Foundation. 

9 http://www.naic.edu/alfa/ 



ing the lower velocity and latitude HI emission and thus 
also do not provide adequate information on the disk/halo 
interface. 

The structure of this paper is organized as follows. We 
start in Section 2 by briefly introducing the consortium 
for Galactic studies with the Arecibo L-band Feed Array 
(GALFA) and emphasizing the importance of the Arecibo 
telescope for large-scale Galactic surveys. Section 3 de- 
scribes the observing strategy adopted for Galactic HI 
spectral-line surveys with ALFA. In particular, the fast 
scanning mode utilizing the basket-weave technique, and 
a newly developed method for obtaining reference HI spec- 
tra by multiple frequency switching, are applicable to ra- 
dio telescopes in general. In Section 4 we summarize major 
previous observational studies of HI in the Galactic anti- 
center. We then present our results in Sections 5 and 6, 
and provide extensive discussion of observational and the- 
oretical issues in Sections 7 and 9. We summarize our 
results and point out future studies in Section 10. 

2. GALFA — BACKGROUND AND FUTURE 

GALFA is a worldwide consortium of scientists (cur- 
rently with 86 members) interested in undertaking large- 
scale Galactic surveys with the Arecibo L-band Feed Ar- 
ray 9 (ALFA). GALFA members have been working to- 
gether with the National Astronomy and Ionosphere Cen- 
ter in planning and executing large surveys to ensure their 
long-reached usage and availability by the greater astro- 
nomical community. The consortium consists of three sub- 
consortia for spectroscopic HI, radio recombination line, 
and continuum surveys. It is anticipated that data prod- 
ucts from GALFA surveys will be archived and easily ac- 
cessible by the astronomical community. 

Galactic surveys with ALFA have greatly improved an- 
gular resolution in comparison with previous large-scale 
single-dish surveys. For example, the all-sky HI survey by 
Bell Laboratories (Stark et al. 1992) had angular resolu- 
tion of 2° , while the HI all-sky survey with the Dwingeloo 
radio telescope (LDS; Hartmann & Burton 1997) had an- 
gular resolution of 0.6°. GALFA's angular resolution of 
only 3.5' will be especially appreciated at high Galactic lat- 
itudes where high resolution, interferometric observations 
are still rare. At lower latitudes several recent, high reso- 
lution, surveys of the Galactic plane exist, or are underway 
(Taylor et al. 2003; McClurc-Griffiths et al. 2005), with a 
typical angular resolution of 1-2' and a brightness sensi- 
tivity of 1.5-3 K (per 0.8 km s _1 wide velocity channels). 
To provide complete information about the large-scale HI 
structure these surveys combine data from interferometer 
and single dish telescopes. 

The Arecibo telescope has several unique advantages 
over all other radio telescopes: (i) its unmatched collecting 
area, together with a large gain and low-noise receivers, re- 
sults in an excellent and unrivaled filled surface-brightness 
sensitivity; (ii) this is the only radio telescope that can 
simultaneously provide both an high angular resolution, 
usually provided exclusively with interferometers, and a 
continuous range of spatial scales (including the largest 
scales, or so-called zero spacings), usually provided only 

Ionosphere Center, operated by Cornell University under a cooperative 
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with single dish telescopes. This unique combination of 
high surface-brightness sensitivity and accurate mapping 
of all angular scales above 3.5' makes the Arecibo tele- 
scope extremely valuable for Galactic science. In addition, 
ALFA, a cluster of seven receivers, has increased Arecibo's 
surveying (mapping) speed by a factor of seven, making 
the telescope a vastly more efficient system for covering 
large areas than an interferometer/single dish combination 
(Dickey 2002). In the fastest mapping mode commonly 
used for GALFA projects, the mapping speed is ~ 9.4 
square degrees per hour. This results in a rms bright- 
ness sensitivity per beam of 0.1 K (per 0.72 km s _1 wide 
channels and after averaging both polarizations). For a 
spectral feature with a FWHM of 10-20 km s _1 our min- 
imum (3-er) detectable HI column density is (5-11) xlO 18 
cm~ 2 . 

This paper focuses on only one of major science drivers 
for GALFA — the Galactic disk/halo interface region. The 
whole science case for Galactic HI spectral-linc surveys 
is rich and diverse, full description can be found in the 
GALFA White Paper 10 . In addition, there are several 
projects planned for continuum and radio recombination 
line surveys which will be presented elsewhere. Currently, 
several different scientific HI projects are underway, some 
focusing on relatively limited areas of the sky (for study- 
ing, for example, particular molecular clouds and their 
possible connection with atomic gas) , while some others re- 
quiring the coverage of the whole sky visible from Arecibo 
(the study of halo cloud properties, for example). The 
general survey philosophy has been to start with smaller, 
scientifically interesting regions, apply exactly the same 
observing parameters (velocity resolution, scanning mode 
etc.), and eventually combine all observations to build the 
whole "jigsaw puzzle" of the Galactic HI sky observable 
from Arecibo (declination range —1° to 38°). 

While large-scale Galactic HI surveys with ALFA will 
address numerous outstanding questions regarding Galac- 
tic science, their importance for other scientific areas is 
also large. For example, X-ray observations need to be 
corrected for Galactic absorption using HI column den- 
sity measurements. The small-scale structure (< 30') in 
the neutral Galactic layer can alter significantly interpre- 
tation of X-ray data and needs to be accounted for (Breg- 
man et al. 2003; Lockman 2004). Also, for the interpreta- 
tion of UV absorption spectra, HI emission spectra with 
as high angular resolution as possible, approximating a 
pencil-beam toward a background source, are necessary 
(Wakker et al. 2003). 

The large-scale HI Galactic survey with ALFA started 
in early 2005. Prior to this, in September-November 2004, 
short precursor observations were undertaken to test the 
instrumentation, observing techniques, and data reduction 
tools. One of those test runs focused on a region of inter- 
est for studies of the Galactic disk/halo interface region, in 
the Galactic anti-center, well known to harbor a wealth of 
high velocity gas. As a demonstration of the quality of HI 
images obtained with ALFA we present here observations 
of this region and preliminary results. 

10 http:/ /www. naic.edu/alfa/galfa/galfa_docs.shtml 

11 http://www.naic.edu/alfa/memos/ 

12 http://www.naic.edu/alfa/memos/ 



3. GALFA OBSERVATIONS AND DATA PROCESSING 

The observations presented in this paper were under- 
taken in October and November of 2004, as a precursor 
observing run of the GALFA Consortium with the newly 
installed ALFA system. Our main aim was testing of 
the newly developed spectrometer (GALSPECT) and the 
proposed observing strategy. ALFA was installed at the 
Arecibo telescope in April 2004. 

ALFA is a cluster of seven state-of-the-art dual polariza- 
tion receivers operating at the frequency range 1225-1525 
MHz. ALFA was constructed at the Australia Telescope 
National Facility. Each of seven beams has a slightly el- 
liptical shape, with the observed FWHM being relatively 
constant at 3.13' x 3.57' (in Az and za) (Heiles 2004a 11 ; 
Deshpande et al. in preparation). The on-axis gain is 10.7 
K Jy _1 for the central beam and 8.3, 8.5, 7.8, 8.3, 8.7 and 
8.3 K Jy^ 1 for the outer six beams, respectively (Heiles 
2004b 12 ). The beam efficiency (with the first sidelobe) 
ranges from 0.80 for the central beam, to 0.69-0.77 for the 
outer beams. The observed beam properties agree very 
well with the theoretical predictions by Cortes-Mcdellin 
(2002). The central beam has the smallest and the most 
circularly symmetric sidelobe; the outer beams have asym- 
metric sidelobes. 

3.1. GALSPECT 

GALSPECT is a dedicated spectrometer for Galactic 
HI surveys. This newly developed FPGA-based system 
was specifically designed for GALEA'S needs. GALSPECT 
takes simultaneously two spectra: a science spectrum cov- 
ering 7.14 MHz with 8192 frequency channels, resulting 
in the fixed velocity resolution of 0.18 km s -1 , and a cal- 
ibration spectrum covering 100 MHz with 512 frequency 
channels. The usable velocity range after data processing 
is -700 to 700 km s" 1 . 

GALSPECT performs the frequency analysis by means 
of a polyphase filter bank (Crochiere & Rabiner 1983). 
This is an extension of the segmented-FFT technique 
(Romncy 1995) which improves greatly the channel sep- 
aration (down to —80 dB) at a small increase in com- 
putational complexity (less than a factor of two). This 
means that there is essentially no 'spectral leakage', which 
results in superb interference isolation. The data are dig- 
itized using 8-bit A/D converters. With this many bits, 
the clipping correction is small and we do not apply it. 

3.2. Observing strategy 

The adopted observing strategy takes advantage of the 
basket-weave (or meridian-nodding) scanning. We observe 
always at the meridian by driving the telescope only in 
zenith angle. Each day a single basket-weave scan is com- 
pleted. As the Earth rotates a scan draws a zigzag pattern 
in the RA-Dec coordinate frame. On consecutive days we 
obtain adjacent scans, shifted by ~ 1.79', covering even- 
tually the whole region of interest. To cover the region 
about 30° x 5° in size that is presented in this paper, we 
required 2 hours x 8 days. The basket- weave scanning 
technique has been commonly used in the past for all-sky 
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Fig. 1. — An image of the Galactic anti-center region from the Leiden- Dwingello survey (Hartmann & Burton 1997), at the LSR velocity of 
—22.67 km s -1 , emphasizing the filamentary structure seen at intermediate velocities. The white box outlines the region covered in GALFA's 
precursor observations. 



continuum surveys (for example see Haslam et al. 1982). 
It has several important advantages when imaging large 
areas of the sky: fast coverage of a large area of the sky, 
the beams always have the same orientation, inter-woven 
scans have many crossing points that allow a fine gain 
adjustment, and transit observations limit contamination 
from position-dependent sidelobes (which at Arecibo can 
be significant). This very fast observing mode scans the 
sky at a rate ~ 5.7 times faster than the sidereal rate, but 
covers every piece of the sky twice. The resulting integra- 
tion time is 2.45 sec for each Nyquist pixel (1.8' x 1.8' in 
size). The separation between ALFA beams is ~ 1.79' on 
the sky, close to Nyquist sampling, and it does not depend 
on the scanning speed. Practical details are given in the 
"GALFA HI User's Guide", Krco et al. (2006) 13 . 

Obtaining reference spectra is one large challenge 
GALFA HI surveys have to deal with. While mapping 
large-scale Galactic structure it is impossible to find an 
'empty' point on the sky that will serve as a reference posi- 
tion. Commonly, the frequency switching technique is used 
when performing large-scale Galactic surveys. However, at 
many single-dish radio telescopes, and especially at the 
Arecibo telescope, the classic frequency switching tech- 
nique introduces baseline problems because of frequency 
structure in the baseline produced by reflections. We have 
developed a new technique for obtaining a reference spec- 
trum, called "The Least-Squares Frequency Switching" 
(LSFS, Hcilcs 2005a 14 ). This technique derives the IF 
bandpass correction from a set of observations obtained 
with 7 LO frequencies. At the beginning and end of each 
day's basket-weave scan we spend about 10 minutes to 
obtain LSFS data. The LSFS technique is based on cy- 
cling among 7 different LO frequencies and using the least- 
squares technique to separate the IF gain from the RF 

13 http:/ /www. naic.edu/alfa/galfa/galfa_docs.shtml 

14 http://www.naic.cdu/alfa/galfa/galfa_docs.shtml#obs_data 

15 http://www.naic.edu/alfa/memos/ 



spectrum. LSFS is performed separately on the 100-MHz 
calibration spectrum and on the 7.14-MHz science spec- 
trum. Broad frequency structure in the IF bandpass and 
baseline of the science spectrum are derived from those 
portions of the calibration spectrum that lie outside the 
science spectrum's band, and appropriate corrections are 
applied to the science spectrum. 

3.3. Data reduction 

A special purpose suite of IDL programs has been de- 
veloped for the reduction of GALFA HI data, primarily 
by C. Heiles and J. Peek. The data reduction consists of 
five stages: bandpass correction, temperature calibration, 
cross-point correlation, fixed pattern noise reduction, and 
gridding. During the first stage, spectra from all beams 
arc corrected for the IF gain, using the baseline as derived 
from the wideband calibration spectrum. We also correct 
here for a fast ripple that arises from reflections in the 
optical fiber that connects the control room with the feed 
cabin. 

In the second stage of data reduction, spectra from all 
beams are calibrated into antenna temperature units using 
a fiducial conversion factor, which is later greatly refined. 
After this, crossing points of all beams are found in the 
data (this is typically a large number ~ 420 / cos(Dec) per 
square degree) , crossing spectra are compared and the rel- 
ative point-to-point gain corrections are determined. At 
the fourth stage, baselines of the beams are partially cor- 
rected for the fixed pattern noise (Heiles 2005 15 ) using a 
beam-averaging process. 

Finally, the last stage of data reduction includes grid- 
ding of calibrated and gain-corrected spectra from all scans 
into a spectral-line data cube. We use a program that as- 
sumes a Gaussian convolving kernel. Prior to gridding all 
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Fig. 2. — An example of interesting structure found in the anti-center data cube: cloudy structure at b = 3°. A chimney-like extension is 
noticeable around I ~ 178° which may be connected with numerous small clouds at the velocity of ~ —40 km s _1 . To enhance weak features 
all pixels with Tb > 30 K have been masked out. 
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Fig. 3. — An example of interesting structure found in the anti-center data cube: three remarkable small HI clouds at b = 18°, I = 186-188°, 
at the LSR velocity —21 km s _1 . All pixels with Tb > 12 K have been masked out. 
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Fig. 4. — An HI velocity profile for the low-velocity cloud at (l,b) =(187°, 18°) as obtained with GALFA and from the LDS data. This is 
the middle cloud seen in Figure 3. These profiles clearly demonstrate that the visibility of emission from a small feature at a velocity of —20 
km s _1 is dramatically improved by a factor of 100 smaller beam solid area achieved by GALFA. 
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data are organized in a sparse matrix form which provides 
a memory-effective way of addressing a huge amount of 
data in IDL. The final post-gridding angular resolution 
is 3.5'. In this particular data cube, a 1-a noise level is 
0.16 K, per 0.74 km s _1 wide velocity channel. The final 
data cube was at the end absolute-calibrated to match the 
brightness temperature scale of the LDS data for this re- 
gion. To achieve this we scaled all spectra by a constant 
factor of 0.84. 

3.4. Effects of stray radiation 

The stray radiation properties of the Arecibo telescope 
are not well measured and understood. Above, we gave 
the main-beam plus first-sidelobe efficiencies; these vary 
from ~ 0.7 to 0.8, which leaves a substantial residual beam 
response from sidelobes. A major question is: do these 
sidelobes reside at large angles from the main beam - 
normally referred to as "stray radiation" — as they do for 
smaller telescopes having significant blockage from struc- 
tures like feed legs, or are they more concentrated toward 
the viewing direction? 

We have performed some rough tests to examine this 
question. We compared emission profiles obtained in the 
Millennium survey (Heiles & Troland 2003a) with the LDS 
data, which are corrected for most effects of stray radia- 
tion. For data taken at modest zenith angles, below ~ 17°, 
the Arecibo and LDS profiles have similar shapes in the 
line wings, where the effects of stray radiation are most 
noticeable. Moreover, our comparison of Millennium sur- 
vey profiles of the same position taken at different times 
of year shows excellent agreement. The absence of signif- 
icant response at large angles from beam center can be 
understood physically: most of the blockage at Arecibo is 
produced by structures that are much larger than a wave- 
length, so their diffraction patterns tend to be concen- 
trated in a fairly tight angle around the beam center. 

However, for pointing directions having zenith angle 
> 17°, the agreement between LDS and Arecibo profiles is 
not so good. Thus, for such directions, stray radiation at 
Arecibo does contribute artifacts in the line wings. This 
division line near zenith angle 17° makes sense because 
above that zenith angle the feed illumination pattern falls 
beyond the edge of the spherical primary reflector and onto 
the surface shield placed around its rim (so called skirt), 
which by its very nature produces far-out sidelobes with 
the consequent stray radiation. A project for the future is 
to pursue the measurement of this stray radiation in more 
detail. 

Given this position-dependent stray radiation, we must 
ask how significant it is for our desired science. The at- 
tractive feature of mapping HI with the Arecibo telescope 
is the good angular resolution, compared to the LDS, com- 
bined with the high surface brightness sensitivity. Stray 
radiation has essentially no scientific impact on this at- 
tractive feature. The telescope's stray radiation response 
comes from weak sidelobes that cover large angles in the 
sky, and these cannot produce structure on angular scales 
measured below a few degrees. In short, changes in the 
measured HI profile on angular scales of a few degrees can- 
not possibly be produced by the convolution of these weak, 
broad sidelobes with structure on the sky because even for 
a point-source sky feature the convolution produces a re- 



sponse that has the angular scale of the sidelobes. 

We conclude that on scales below of a few degrees, angu- 
lar structure — i.e., change in the HI profile — is not affected 
by distant sidelobes. However, on the same scales the over- 
all shape of our profiles are sometimes susceptible to con- 
tamination by stray radiation. These overall shapes need 
to be corrected by referring to a stray-radiation-corrected 
HI survey, of which the best is the Leiden- Argentine-Bonn 
(LAB) survey (Kalberla et al. 2005a; Hartmann & Bur- 
ton 1997; Bajaja et al. 2005; Arnal et al. 2000). This 
survey has ~ 36' resolution, and when our maps cover a 
reasonable area of a few square degrees we can correct our 
profiles by using this survey as a standard. We have not 
yet implemented such a correction procedure. 

4. HISTORY OF STUDIES OF ANOMALOUS GAS IN THE 
ANTI-CENTER REGION 

As mentioned in Section 2, the first GALFA observa- 
tions focused on a region in the Galactic anti-center. The 
Galactic anti-center is well known as a region with a large 
amount of gas at forbidden velocities and hence offers 
the perfect case-study for Galactic disk/halo interfaces. 
Weaver (1970) was the first to point out interesting high- 
velocity, jet-like features in the Galactic anti-center region 
which were apparent at some longitudes and appeared to 
be connected to the low-velocity gas in the Galaxy. Sev- 
eral authors have studies these features further, for the full 
list of references please see Tamanaha (1997). 

The most prominent structure in the anti-center region 
is the Anti-center Shell (ACS), discovered by Heiles (1984), 
and studied by several authors. This large supershell, with 
an angular diameter of almost 45°x25°, is moving very 
fast (velocity > 70 km s _1 ). However it appears to have 
stopped expanding. Several authors, including Kulkarni 
et al. (1985), have favored the explanation that the ACS 
was formed by an impact of a chain of high velocity clouds 
onto the Galactic disk. Tamanaha (1994) pointed out that 
the ACS has a wealth of filamentary structure at interme- 
diate velocities, —25 to —90 km s _1 . By using a special 
filtering technique he was able to isolate seven filaments 
and propose a model in which filaments are located on the 
surface of a sphere with a radius of about 25°, centered on 
I = 185°, b = 0°. These filaments are large, having a typi- 
cal length of 10°-40°, and appear at mostly constant lon- 
gitudes, and have typically FWHM of 10-20 km s _1 . The 
chain of HVCs in the anti-center region, often referred to 
as clouds ACI, ACII etc., has been known for a long time. 
It is still unclear though whether these clouds are related 
to the ACS. 

5. THE ANTI-CENTER FILAMENTS AND 'COLD' CLOUDS 

This 'pilot' HI spectral-line data cube is centered on 
the Galactic anti-center region, I ~ 180°, and covers a 
range of latitudes, b =0°-23°, with Vi SI = 90 to -200 
km s . The original velocity resolution was 0.18 km s _1 , 
we have smoothed spectra by four channels resulting in 
the final velocity resolution of 0.74 km s _1 . In the veloc- 
ity domain the data cube covers several distinctly differ- 
ent environments: Galactic HI, intermediate velocity gas 
at Vi sr ~ —50 km s _1 , a portion of the Outer Arm HVC 
complex at Vi SI ~ — 80 km s _1 , and a typical compact 
high velocity cloud (CHVC) at Vi sr 115 km s" 1 . This 
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Fig. 5. — (left) An image obtained with ALFA showing filamentary structure at low forbidden velocities, around —22 km s _1 . The grey-scale 
range is 1 to 21 K with a square-root transfer function to enhance weak features, (right) A zoom in on the weak filamentary structure at 
b ~ 18° with a stretched intensity scale. 



CHVC was cataloged as WvW215 by Wakker & van Woer- 
den (1991) and as CHVC186+19-114 by Braun & Burton 
(1999). 

Figure 1 shows a large view of the region, obtained from 
LDS. The white box included in this figure encloses the 
area covered in GALFA's precursor observations. Our data 
offer one of the very first views of the high-latitude HI dis- 
tribution at the high angular resolution of ~ 3.5' with 
high surface brightness sensitivity. This view reveals sev- 
eral new and interesting phenomena, which we introduce 
in the following subsections. The compact HVC is intro- 
duced in Section 6. 

5.1. Low-velocity clouds or LVCs 

The region contains numerous discrete HI clouds at neg- 
ative velocities of V\ SI ~ —20 to —40 km s _1 , shown in Fig- 
ures 2-6. Throughout this paper we will call these clouds 
'low- velocity clouds' (LVCs). Although LVCs are in the 
forbidden velocity range as the expected LSR velocity for 
gas in the Galactic anti-center and under Galactic rotation 
is ~ km s _1 , velocities up to 20-30 km s _1 are typically 
allowed to account for turbulent motions (Wakker 2001). 
Most of LVCs have |Vi sr slightly smaller than, or near the 
low velocity end of, what is traditionally assumed for inter- 
mediate velocity clouds (e.g. IVCs are defined as having 
|M S r| > 30-40 km s" 1 ; Wakker 2001, Richter et al. 2003). 
Frequently, LVCs appear to kinematically follow the main 
Galactic disk structure. 

This type of HI cloud is particularly common at the 
lower Galactic latitudes in our region. An example at 
b = 3° is shown in Figure 2, where numerous, discrete, 
narrow-lincwidth HI clouds are easily noticeable. Some 
of these discrete clouds are connected in velocity to lower 
velocity gas; the strongest example is at £ = 178°. In 
this particular data cube HI clouds like these are not com- 
mon at intermediate Galactic latitudes, b =~ 7 — 15°, 
while they become again very frequent at high latitudes, 
b > 15°. Clouds at high Galactic latitudes are especially 
easy to notice, and most of them appear to be associated 



with faint, diffuse filamentary structures. 

A spectacular example of LVCs at b = 18° is shown in 
Figure 3: three small HI clouds (the smallest cloud has an 
angular size of only 4.5'), with extremely narrow velocity 
lincwidths (~ 2 km s _1 ) are noticeable at distinctly differ- 
ent LSR velocities from those of the bulk HI disk emission. 
However, the clouds clearly follow the distribution of the 
disk HI with longitude, and are probably affected by the 
general Galactic rotation. Figure 4 shows spectra of one 
of the three clouds as obtained with GALFA (left panel) 
and from LDS (right panel). The LDS profile only barely 
hints at the —20 km s _1 feature, which is so prominent in 
the GALFA profile. Clearly these HI clouds are too small 
to be visible in low-resolution surveys and were missed. 

5.2. Low-velocity filamentary structures 

The region has prominent filamentary structures, espe- 
cially near Vi sr = —25 km s _1 — covering roughly the 
same velocity range as that of the discrete clouds. These 
filaments are clear in Figure 5 (left), which shows a veloc- 
ity channel at Vi sr = —22.1 km s~ x . The filaments may be 
related to the larger filaments noted previously in the anti- 
center region and studied extensively by Tamanaha (1994, 
1997); however the character of these filaments changes 
dramatically when viewed with our full angular resolution. 
Figure 5 (right) zooms in on the black square in the left 
panel. With this higher magnification, these filaments are 
peppered with numerous embedded HI clouds. Most of 
the clouds are small, ~ 5'-10' in size, and very difficult to 
notice in previous low-resolution HI images. 

5.3. Cloud morphology 

Several LVCs have especially interesting morphology 
suggestive of core/envelope structure which is indicative of 
the multi-phase medium. Figure 6 (left) shows an example 
of a possible connection between the disk HI and HI clouds 
at Vi sr = —38 and —50 km s _1 — two weak HI 'bridges' 
(Tb ~ 1 K) are noticeable extending from the disk to the 



8 



S. Stanimirovic et al. 




Fig. 6. — Two examples of interesting cloud morphologies found in the anti-center data cube. To enhance weak features all pixels with 
Tb > 30 K have been masked out. (left) An example of a larger cloud at (l,b) = 186.9,16.4° at a velocity of —38 km s _1 which shows 
weak elongations (tails), with Tb ~ 1 K, at both lower and higher velocities. Similarly, the intermediate-velocity cloud at a velocity of —50 
km s has a tail (Tb ~ 1—1.5 K) connecting it with the disk HI emission. Contours arc at 1.0, 1.25, 1.5 and 3.0 K. The grey-scale range is 
—2 to 12 K with a linear transfer function, (right) A cloud at (I, b) = 183.8, 10.4°, at a velocity of —43 km s _1 which appears to have a weak 
tail trailing all the way to —80 km s _1 . Contours range from 5 (0.75 K) to 10-ct (1.5 K). 



intermediate velocity gas, suggesting a possible physical 
connection. Both clouds and bridges could also repre- 
sent segments of a larger, shell-like structure. There is 
no evidence, however, that the size of the shell- like feature 
changes systematically with latitude, as it would in the 
case of an expanding shell. The cloud 183.81+10.37, indi- 
cated in Figure 6 (right) with the black arrow, is another 
interesting example — it has a faint 'tail' (Tb ~ 1 K) that 
trails toward the intermediate-velocity gas at Visr ^ 

-80 

km s _1 . 

5.4. Properties of low-velocity clouds 

To avoid confusion with Galactic emission we have se- 
lected a small set of 12 clouds with b > 10° which stand 
out as obvious, well-separated features. These clouds are 
distinct, well separated from the Galactic plane emission, 
and have negative deviation velocities —42 < Vd ov < —25 
km s _1 . Wakker (2001) defined the deviation velocity as 
the minimum value of the deviation of the cloud's radial 
velocity from Galactic rotation. In this framework, clouds 
with negative Vj sr and negative Vdev appear to be moving 
toward us too fast relative to what is allowed by Galactic 
rotation. This set of clouds is by no means complete. Our 
aim is to demonstrate the existence of small, cold HI clouds 
in the outer Galaxy, at forbidden velocities and most likely 
located in the Galactic disk/halo interface region, as well 
as to describe their basic properties. With a larger amount 
of collected data in the near future we will be able to do 
a statistically more meaningful analysis. 

Figure 7 shows velocity profiles through the centers of 
selected HI clouds. HI clouds/condensations easily stand 
out as distinct peaks with a central velocity in the range 
~ —25 to ~ —42 km s _1 . To measure the clouds' basic 
properties we have fitted each velocity profile with a few 
Gaussian functions. Results from this fitting procedure 
are summarized in Table 1 where we provide the range of: 
cloud angular sizes, Vj sr , velocity FWHM (Av), the peak 
brightness temperature (TB.max), and the peak HI column 
density (iV(HI) poa k). In Figure 8 central positions of se- 
lected clouds are shown with crosses on a I — V\ ST image. 
Each pixel in this image represents the peak brightness 
temperature over the latitude range b =15°-19°. This im- 
age depicts almost continuous extent of cloudy structure 
between the disk HI emission and the intermediate velocity 



gas at ~ —50 km s _1 . Many of LVCs appear to kinemat- 
ically follow the disk HI, but 2-3 clouds with the most 
negative velocities could be related to the intermediate- 
velocity gas. Several small clouds at positive LSR veloci- 
ties, around +20 km s _1 , are also noticeable in this figure. 

Most of LVCs must be very cold. The velocity FWHM 
for all 12 clouds is in the range 3.0-7.6 km s _1 . The mean 
FWHM is 4.4 km s _1 , and the mean Tk jmax (the kinetic 
temperature in the case there is no nonthermal broaden- 
ing) is only 470 K. 80% of clouds in this (small) sample 
have T k < 450 K and FWHM of 3-4 km s" 1 . The clouds' 
mean central LSR velocity is —33 km s _1 and the mean 
peak HI column density is 2 x 10 19 cm~ 2 . The cloud an- 
gular size ranges from 6' to 12'. Here we focus only on 
well-isolated and compact clouds; but there are also clouds 
with diameters of up to 30' in the data cube which we will 
discuss in a future paper. 

5.5. Properties of low-velocity filaments 

Most of these clouds appear related to diffuse filamen- 
tary features (Figure 4) which are seen as broad line wings 
in the Galactic HI profile. From the fitting of velocity pro- 
files as well as visual inspection, the filaments have a large 
FWHM, typically 20-30 km s _1 , and their peak brightness 
temperature is low, Tb = 1-2 K. Their sizes are of order 
of a few degrees, some filaments may extend beyond our 
maps. Typical HI column densities are (4-12) x 10 19 cm~ 2 . 

6. CHVC186+19-114 

Figure 9 shows HI channel maps of CHVC186+19-114 
which is present in the same data cube at higher negative 
velocities. The top ten panels in this figure show primarily 
a more diffuse, or envelope structure, having asymmetric, 
ring-like morphology. The same panels also show an inter- 
esting 'companion' HI cloud, seen for the first time in our 
observations, and located ~ 50 arcmin to the southwest of 
CHVC186+19-114. The more compact (core) component 
of this CHVC is primarily at higher negative velocities, 
typically between —115 and —122 km s -1 , and is visible 
in the remaining panels of Figure 9. 

6.1. Properties of core/envelope structure 

To obtain physical parameters of the HVC's core and en- 
velope separately, we have fitted each velocity profile with 
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Table 1 

Measured properties of HI clouds and filaments at low and high negative velocities. 





Angular size 


v lsr 




FWHM 


Tb , max 


ivrm) peak 




(arcmin) 


(km s 




(km s -1 ) 


(K) 


(10 19 cm" 2 ) 


Low-velocity clouds 


6-12 


-25 to 


-42 


3.0-7.6 


1.6-4.1 


1-3 


Low-velocity filaments 


a few deg 


-20 to 


-30 


20-30 


1-2 


4-12 


CHVC186+19-114 Core 


30 X 40 


-110 to 


-118 


4-7 


~ 1-4 


~ 1-11 


CHVC186+19-114 Envelope 


30 x 65 


-100 to 


-118 


~ 16 


0.5-1.5 


3-13 


CHVC186+19-114 Companion 


8 


-100 to 


-107 


~ 15 


0.8 


6 
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Fig. 7. — Velocity profiles of 12 small, isolated HI clouds. The thin solid line shows the data, the thick solid line is the combined best 
fit, the dashed line shows the Gaussian component corresponding to HI low-velocity clouds, and the dot-dashed line depicts a broad, diffuse 
component, which is almost always needed to produce a good fit. 
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Fig. 8. — An HI peak brightness temperature image over the latitude range b = 15° to 19°. The grey-scale range is to 17 K, with a linear 
transfer function; pixels with Tb > 25 K have been masked out to enhance weaker features. Crosses show central positions of selected clouds 
in our sample. 
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Fig. 9. — HI images of CHVC186+19-114 at the LSR velocities specified in the top-left corner of each panel. The grey-scale intensity range 
is —0.5 to 2 K with a linear transfer function. 
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Right Ascension (J2000) 



Fig. 10. — HI column density image of the broad velocity component shown in grey-scale. HI column density distribution of the narrow 
velocity component is shown with black contours. Contour levels range from 1 X 10 19 to 1.1 X 10 21 cm , with a step of 2 X 10 . 




Fig. 11. — The velocity field of the narrow velocity component (left) and the wide velocity component (right). 
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Fig. 12. — An example of extremely narrow velocity profiles found at several locations in CHVC186+19-114. This profile was obtained at RA 
07 h 18 m 2.5 s , Dec 31°44' 46.5" (J2000) and was fitted with two Gaussian functions: the narrow function has a FWHM of 3.7 km s" 1 (shown 
with the dashed line), and the wide function has a FWHM of 14.7 km s -1 (shown with the dot-dashed line). The total fit is shown with the 
thick solid line. 
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one or two Gaussian functions. For this exercise, only pix- 
els with a brightness temperature T-q > 3a were selected. 
Regions spatially corresponding to the core structure are 
typically well-fitted with a combination of one narrow and 
one wide Gaussian function. Most of the envelope, as well 
as the HI companion cloud, is well-fitted with a single, 
broad Gaussian function. Figure 10 shows HI column den- 
sities of the narrow (contours) and wide (grey-scale) Gaus- 
sian components, while their velocity fields are shown in 
Figure 11. 

An interesting result is that both narrow and wide veloc- 
ity components show a clear velocity gradient from ~ — 104 
km s _1 in the south to ~ —118 km s _1 at the north. 
This suggests that the same physical mechanisms affect 
and de- accelerate both cold and warm gas in this CHVC. 
The central LSR velocity of the companion cloud is —107 
km s _1 , while its typical FWHM lincwidth is ~ 18 km s _1 . 
The wide component has a typical FWHM linewidth of 
~ 15 — 18 km s , while the narrow component has a 
FWHM of ~ 7 km s _1 . Interestingly, there are several 
small regions in this CHVC where the narrow velocity 
component has much smaller linewidth. As an example 
Figure 12 shows a velocity profile where the colder com- 
ponent has a FWHM of only 3.7 km s _1 , indicative of cold 
gas with the kinetic temperature < 300 K. The narrow and 
wide velocity components of the main cloud have compa- 
rable peak iV(HI) of ~ 10 20 cm~ 2 , while the companion 
cloud has a peak A?" (HI) of ~ 5 x 10 19 cm~ 2 . 

As mentioned earlier, the main cloud of CHVC186+19- 
114 was cataloged previously (Wakker & van Woerdcn 
1991; Braun & Burton 1999) and also partially imaged 
at high resolution with the Arccibo telescope (Burton 
et al. 2001) and with the Westerbork Synthesis Radio Tele- 
scope (de Heij et al. 2002). Burton et al. (2001) found 
that CHVC186+19-114 has a core/envelope morphology 
and fitted its velocity profiles with two Gaussian func- 
tions. Based on the measured exponential AT(HI) pro- 
file as a function of cloud's radius, being interpreted as 
due to a spherical exponential distribution of the HI vol- 
ume density distribution, Burton et al. (2001) estimated 
CHVC186+19-114's distance to be 600-800 kpc. Mal- 
oney & Putman (2003), however, showed that exponen- 
tial AT (HI) profiles could be a natural consequence of an 
external source of ionizing photons and cannot be used 
to derive HVC's distances, de Heij et al. (2002) noted 
a velocity offset between the CNM component, traced in 
intcrferometric data, and the WNM component, derived 
from the total power, Arecibo data. They suggested that 
this systematic offset could be explained by some type of 
external perturbations, one possibility being ram-pressure 
interactions with the external medium. 

6.2. Relation of main and companion clouds 

GALFA observations cover a significantly large area 
around CHVC186+19-114, revealing a 'companion' HI 
cloud. The companion cloud has the central LSR velocity 
very similar to that of CHVC186+19-114's envelope and 
lags behind CHVC186+19-114's core by - 10 km s" 1 . In 
addition, velocity profiles (Tb and FWHM) on the south 
side of the HVC envelope are similar to profiles through 
the companion cloud. This is all suggestive of a possible 
physical association between CHVC186+19-114 and the 



companion cloud. One possibility is that this small cloud 
has been stripped from CHVC186+19-114's envelope due 
to ram pressure interactions with the halo medium. This 
cloud could represent an evidence for the break up of HVCs 
into smaller clouds due to ram pressure (or even tidal) in- 
teractions. Another possibility is that both CHVC186+19- 
114 and the companion cloud are embedded in lower col- 
umn density gas and we are just able to detect the high 
column density peaks of this extended 'complex'. 

The companion HI cloud has an angular size of only 
7' x 9' and could be classified as one of the smallest HVCs 
ever detected. Its size is just at the lower end of the 
size distribution of the mini-HVC population discovered 
recently by Hoffman et al. (2004). 

7. OBSERVATIONAL DISCUSSION: LOW- VELOCITY CLOUDS 

7.1. Distances and HI masses of low-velocity clouds 

A major unknown parameter of LVCs is their distance. 
As we have shown in Section 5.4, 80% of LVCs in our 
sample have T k < 400 K and FWHM of only 3-4 km s" 1 . 
The peak AT(HI) for all clouds is (1-3) xlO 19 cm~ 2 . These 
linewidths and column densities are very similar to those 
of ordinary CNM clouds at |6| > 10° studied through HI 
absorption measurements (Hcilcs & Troland 2003b), sug- 
gesting that the anti-center LVCs are likely to be "true" 
CNM clouds. As mentioned earlier, LVCs often appear 
related to larger, warmer filamentary structures, with a 
typical FWHM of 20-30 km s" 1 . 

In the Galaxy, the CNM and WNM can co-exist only 
over a certain, well-defined range of thermal pressures, 
P m i n to P max [this was first emphasized by Field et al. 
(1969) and most recently re-examined by Wolfire et al. 
(2003)]. This co-existence of cold LVCs with warm fil- 
amentary structures can be used to estimate cloud dis- 
tances. As the ISM is highly turbulent, often turbulent 
pressure is more important than the thermal pressure. 
However, as nicely explained in Wolfire et al. (2003), on 
the scales of typical CNM clouds thermal pressure dom- 
inates, and CNM clouds are expected to be embedded 
in a spatially isobaric medium with their surface layers 
having the same thermal pressure as the ambient warm 
medium. Wolfire et al. (2003) calculated P m i n and P max 
in the Galactic midplane for several Galactocentric radii 
(i? g ). Both quantities decrease with R g . For example, at 
R g = 8.5 kpc, P min = 1960 K cm" 3 and P max = 4810 K 
cm~ 3 , while at R g — 18 kpc, P min = 272 K cm~ 3 and 
P max = 1220 K cm" 3 . 

We can now calculate for our LVCs their volume density 
(n) and thermal pressure (P) using: n oc A^(HI) pca k x D^ 1 
and P oc n7k jmax , where Tk jmax = 21.86Af . Both n and 
P depend on the cloud distance D, also these equations 
assume that clouds are spherically symmetric (please note 
that sheets will have higher pressure, while cigars oriented 
towards the observer will have lower pressure). Estimates 
for both n x D and P x D are given in Table 2, together 
with the estimated maximum kinetic temperature and the 
peak HI column density. We can now compare measured 
values for P x D, for each cloud, with the thermal pressure 
allowed for the CNM and WNM co-existence to estimate 
D. If we assume the typical thermal pressure in the Solar 
neighborhood of 3000 K cm -3 , then clouds in our sam- 
ple have the distance range 0.1-1.3 kpc, with the median 
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Derived (median) properties of HI clouds and filaments at low and high negative velocities. 





Tk,max 


iV(HI) peak 


n X D 


nT k>max x D 




(K) 


(10 19 cm" 2 ) 


(cm -3 X kpc) 


(10 3 K cm- 3 x kpc) 


Low-velocity Clouds 


200-1300 


1.0-3.0 


1.3-5.5 


0.3-4 


Low-velocity Filaments 


9000-20000 


4-12 


0.7-1.5 


10-20 


CHVC186+19-114 Core 


350-1000 


2-11 


0.9 


1 


CHVC186+19-114 Envelope 


5600 


2-13 


0.8 


4 


CHVC186+19-114 Companion 


5000 


5 


3.2 
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distance being 0.3 kpc. The lowest thermal pressure that 
allows the CNM and the WNM to co-exist is ~ 270 K 
cm~ 3 (Wolfire et al. 2003), this places an upper limit on 
the distance to our clouds of ~ 3.6 kpc (this value is the 
median for the whole cloud sample). The highest thermal 
pressure that allows the CNM and the WNM to co-exist 
is ~ 5000 K cm~ 3 (for i? g ;> 8.5 kpc as our clouds are 
in the direction of the Galactic anti-center), this places a 
lower limit on the distance to our clouds of ~ 0.2 kpc. 
If the thermal pressure is, however, > 5000 K cm~ 3 then 
cloud distances could be < 0.2 pc and our clouds may 
be very closeby. Such high thermal pressures are not ex- 
pected to be common in the Galaxy for i? g > 8.5 kpc from 
the theoretical point of view. If present they would in- 
dicate the existence of very short-lived, transient clouds. 
The observed morphology of LVCs, in the form of dis- 
crete, isolated clouds with sharp edges and appearing over 
many velocity channels, does not suggest their transient 
character. However, future direct distance measurements 
are essential to confirm cloud distances. We also caution 
that our distances were estimated under the assumption 
of the local pressure equilibrium at cloud surfaces. Sev- 
eral recent observations and numerical simulations have 
challenged the global pressure equilibrium condition for 
the interstellar gas, however still mildly favoring the local 
pressure equilibrium state (for example see Mac Low et al. 
2005). 

Table 2 also gives PxD for low- velocity filaments, which 
appear higher than for the clouds and are most likely due 
to a non-spherical, sheet-like morphology. This is sup- 
ported by the visual appearance of filaments which have a 
far more elongated morphology than the clouds. 

If we assume for now that the LVCs have distances in the 
range 0.2-3 kpc, then their height from the plane is in the 
range 60-900 pc. As the typically assumed scale height 
(half-thickness) of the Galactic thin disk is Zthin = 100 
pc (Belfort & Crovisier 1984; Kulkarni & Heiles 1988; de 
Avillez 2000), most of the LVCs appear located in the tran- 
sition region between the thin disk and the hot Galactic 
halo. With a typical angular radius of 9', LVCs' linear size 
is 0.5-8 pc, while their HI mass ranges from 3 x 10~ 2 to 
7 M©. If the clouds are embedded in the warm medium 
with a temperature of 6000 K, then the mass- loss rate 
due to evaporation (calculated using eq. [47] in McKee & 
Cowie (1977)) is 1.5 x 10 24 gr yr _1 , and hence the cloud 
evaporation timescale is ;> 100 Myr. 

Several previous observational studies have found cold 
clouds at significant heights above the disk in our own 
Galaxy and in several other spiral galaxies. Crawford 
et al. (2002) obtained Nal and Call observations towards 



the opening of the Local Interstellar Chimney and found 
cold neutral clouds with low negative velocities, located 
at 0.3 < \z\ < 2.5 kpc. They suggested that these clouds 
most likely belong to a scattered population of infalling 
diffuse clouds, consistent with a Galactic Fountain model. 
Howk (2005) discussed evidence for cold, CNM-like, gas 
in the disk/halo transition region of several spiral galax- 
ies based on optical imaging of extra-planar dust. They 
suggested that this gas was most likely expelled from the 
thin disk through a quiescent process. Richter et al. (2003) 
found evidence for H2 in dense, mostly neutral phase, and 
closely linked to CNM clouds in IVCs at even higher z 
heights. These results support our distance estimate, how- 
ever the formation mcchanism(s) for cold clouds in the 
disk/halo transition regions remain under debate (see Sec- 
tion 9.1). 

7.2. Is there a link between our low-velocity clouds and 
Lockman's clouds? 

One of the most important results of this study is that 
GALFA observations clearly demonstrate the cloudy struc- 
ture of the Galactic disk/halo interface region in the outer 
Galaxy. The cloudy disk/halo interface, emphasized re- 
cently by Lockman (2002) in the inner Galaxy, is obvi- 
ously not restricted just to the inner Galaxy and could be 
quite extended radially. We compare here properties of the 
disk/halo interface clouds in the inner and outer Galaxy. 

Recent HI observations with the Green Bank telescope 
showed that the halo in the inner Galaxy is not smooth 
but populated with discrete HI clouds (Lockman 2002). 
These clouds, often referred to as "Lockman's clouds", are 
located ~ 900 pc below the plane. However, they appear 
to kinematically follow the HI distribution in the disk. We 
summarize typical properties of this cloud population in 
Table 3. Frequently, clouds are organized into larger struc- 
tures and connected by diffuse envelopes and/or filaments 
which can have sizes up to ~ 1 kpc. In addition, high res- 
olution observations (Lockman & Pidopryhora 2004) show 
that most Lockman's clouds have significant sub-structure, 
suggestive of the presence of a two-phase medium. 

Stil et al. (2005) discovered an additional set of 17 iso- 
lated HI clouds in the inner Galaxy, with velocities up to 
60 km s _1 greater than the maximum velocity allowed by 
Galactic rotation (see Table 3). Despite their forbidden 
velocities, these clouds follow the main disk emission and 
therefore cannot be classified as high- velocity gas. Located 
within the Galactic thin disk, at \z\ <J 80 pc, these fast- 
moving clouds demonstrate that interstellar clouds can 
have very large random velocities along the line-of-sight. 
Stil et al. (2005) suggested that the fast-moving clouds, at 
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Table 3 

Median properties of 'different' types of HI clouds in the Galactic disk/halo interface region. 



Reference Number of I b V dev Size FWHM AT(HI) peak D \z\ M m 

clouds (°) (°) (kms- 1 ) (pc) (kms-') (xl0 1H c m - 2 ) (kpc) (pc) (Mg) 

Lockman (2002) 38 21T29 -4 to -12 15 to 40 a 24 12 2 O 950 < 50 

Stil et al. (2005) 17 18-67 -1.3 to 1.3 25 to 60 ~ 10 b 6 20 8 80 60 

This work 12 183-188 12 to 23 -25 to -40 0.5-8 4 2 0.2-3 60-900 0.03 - 7 



"One half of clouds in this study has Vd ev < 15 km s 1 . Median Vd, 
6 Most of the clouds in this study appear significantly elongated. 



\z\ <J 80 pc, and Lockman's clouds, at \z\ ~ 950 pc, may 
belong to the same population of clouds that is widely 
spread throughout the Galaxy. Both works focused exclu- 
sively on clouds at tangent points in the inner Galaxy and 
therefore have well-determined cloud distances. 

Table 3 lists properties of LVCs in the outer Galaxy 
to facilitate their comparison with HI clouds in the inner 
Galaxy. This table gives: the number of clouds in each 
study, the range of Galactic longitude and latitude, cloud 
deviation velocities (as defined in Section 5.4), median lin- 
ear sizes, FWHM, peak HI column densities, typical dis- 
tances, heights above the plane, and estimated HI masses. 

Several potentially important differences are noticeable. 

(1) It is striking that our LVCs, which lie outside the So- 
lar circle, are on average smaller, colder, and less massive 
than clouds in the inner Galaxy. For our selected sample, 
most clouds have FWHM of 3 — 4 km s" 1 , a factor of three 
smaller than that of Lockman's clouds. Angular sizes of 
LVCs are in the range 6'-12' (with a few exceptional ex- 
amples having a size of only 4'), while Lockman's clouds 
have a median size of <~ 11'. 

Linear sizes and HI masses of LVCs are of course uncer- 
tain as they depend on cloud distances. However distances 
of 6-8 kpc, which would bring HI masses close to those of 
Lockman's clouds, would imply the thermal pressure of 100 
K cm" 3 which is below the minimum pressure at which the 
CNM can exist in thermal equilibrium (Wolfire et al. 2003). 
Cloud angular sizes in our study though could be partially 
a selection effect, as we have focussed on well-isolated com- 
pact clouds for which Lockman's angular resolution is less 
adequate. 

(2) | Vdev | of LVCs is on average higher than that of Lock- 
man's clouds. The median deviation velocity for Lock- 
man's sample is 13 km s" 1 , while LVCs have (|Vdev|) = 30 
km s" 1 . If real, this difference could be important, sug- 
gesting that interface clouds at larger galactocentric radii 
(> 8 kpc) have larger deviation velocities than clouds at 
smaller radii (<J 5 kpc). Again, the difference in |Vdev| is 
at least partially due to our selection criterion — to avoid 
confusion with the disk emission we have selected clouds 
with | Vdev | <; 25 km s" 1 , but clouds with lower |Vd ev | 
also exist in the data set. Larger cloud samples from the 
GALFA surveys will be able to clarify whether clouds with 
larger Vdev are more common in the outer Galaxy. 

(3) It is also important to notice that LVCs have nega- 
tive Vdev , meaning they are moving toward us at velocities 
too fast to be explained by the Galactic rotation. Lock- 
man's clouds, however, have systematically positive Vdev, 



= 13 km s" 1 . 



meaning they are moving away from us with velocities too 
fast to be explained by the Galactic rotation. 

In conclusion, if our clouds would have distances of 6— 
8 kpc then their properties would be similar to those of 
Lockman's clouds, but this would imply the thermal pres- 
sure below what is traditionally allowed for the CNM. If, 
however, our clouds are at distances <J 2 — 3 kpc, then 
they are colder, smaller, and less massive than Lockman's 
clouds. Taking into account various selection effects and 
different resolution provided by the Arecibo and Green 
Bank telescopes, clouds in the inner and outer Galaxy may 
belong to the same population with a range of properties 
(z, sizes, HI masses). However, distinctly different forma- 
tion mechanisms can not be excluded either at this stage. 
Halo clouds with sizes similar to those of our LVCs were 
reported recently by Kalberla et al. (2005b) and Dedes et 
al. (in preparation). Using the Effclsberg radio telescope 
these authors found about 20 HI clouds at R g — 13 — 17 
kpc and z = 3 — 5 kpc. However, when viewed at high 
resolution with the Very Large Array, these clouds resolve 
into smaller cores with a FWHM of 3-7 km s" 1 , a peak 
HI column density of a few xlO 19 cm" 2 , and diameters of 
a few pc. 

8. OBSERVATIONAL DISCUSSION: ULTRA-COMPACT HVCS 

The CHVC studied in this paper, CHVC186+19-114, is 
fairly typical of the ^250 currently cataloged CHVCs with 
diameters generally on the order of 20'-80' (Putman et al. 
2002; de Heij et al. 2002). The resolution and sensitiv- 
ity of GALFA observations offer an unique opportunity to 
detect a large number of newly identified ultra-compact 
HVCs (UCHVCs), or HVCs with diameters £ 20'. The 
number of these UCHVCs that have been identified to 
date remains limited. The smallest was found by Brians 
& Westmeier (2004) and has an angular size of 4.5' and a 
peak HI column density of 2 x 10 20 cm" 2 (Brims & West- 
meier 2004). Slightly larger UCHVCs (> 9'; referred to 
as mini-HVCs) were found by Hoffman et al. (2004) with 
much lower peak column densities (< 10 19 cm" 2 ); a simi- 
lar UCHVC was mapped by Richter et al. (2005). These 
UCHVCs have been completely missed in the lower resolu- 
tion surveys and it remains unclear if they will follow the 
power law distribution of flux density, column density and 
size found for other HVCs, or if they represent a distinct 
population of objects. 

The companion cloud in Figure 8 has a size of only 7' x 9' 
and a linewidth of 18 km s" 1 and can thus be classified 
as one of the smallest HVCs known. These UCHVCs may 
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represent a variety of objects. They may be the break- 
ing up of a large cloud as it is stripped via a combination 
of tidal and ram pressure forces. The UCHVCs may also 
represent high column density peaks in a mostly ionized or 
lower column density high- velocity complex. In this case 
the UCHVC may represent the beginnings of a halo cloud 
cooling out of the hot diffuse medium. 

The close match between the position and velocity 
of the UCHVC and the CHVC in our data does indi- 
cate they are directly related. This is similar to the 
CHVCs in the vicinity of the Magellanic Stream that 
show a spatial and kinematic link to the main structure 
of the Stream (Putman et al. 2003; Brims et al. 2005). 
Also, properties of the companion cloud (angular size and 
FWHM) arc similar to those of HI condensations (cores) 
observed with the Westerbork Synthesis Radio Telescope 
in the case of CHVC120 - 20 - 443 (de Heij et al. 2002). 
CHVC120 - 20 - 443 is one of the only two known CHVCs 
with exceptionally broad lincwidths in its CNM (core) con- 
densations, and de Heij et al. (2002) suggest that it may 
be undergoing ram-pressure or tidal stripping. Finally it 
remains possible that the UCHVCs represent small dark 
matter halos as proposed for the CHVCs (Braun & Burton 
1999). Future GALFA observations of larger areas on the 
sky will be important for establishing the spatial distribu- 
tion and properties of UCHVCs. 

9. DISCUSSION: THEORETICAL ISSUES 
9.1. Galactic Fountain and ballistic clouds 

One of the most common models for the interaction be- 
tween the Galactic disk and the halo is based on Galactic 
fountain flows (Shapiro & Field 1976; Houck & Brcgman 
1990). In this model, hot gas from the disk is pushed up 
by the stellar activity in the form of buoyant outflows, it 
travels through the halo, cools down, and as a result of 
thermal instabilities rains back onto the disk in the form 
of cold clouds. In their best model for low-temperature 
fountains (model 3; for a halo temperature of 3 x 10 5 
K and a volume density 10~ 3 cm~ 3 ) Houck & Brcgman 
(1990) suggested that about 70 cold halo clouds should be 
present above a disk area of 400 pc 2 , having a typical size 
of <~ 3 pc, a temperature of 500 K, and a column density of 
2.5 x 10 18 cm -2 . A peak in cloud distribution is expected 
at the location of cloud formation, z ~ 1-2 kpc, however 
if drag is important then the density of clouds at z < 0.5 
kpc could be significant too. The expected radial veloci- 
ties for fountain clouds are in the range from —100 sin & to 
+50 sin b km s _1 . 3-D simulations of a Galactic fountain 
by de Avillez (2000) also show the formation of numerous 
small cloudlets, with sizes down to their resolution limit of 
1.25 pc, resulting from a breakup of collimated structures 
and also larger clouds being swept up by blast waves. 

Typical sizes of LVCs are surprisingly similar to what is 
expected for fountain clouds. In addition, at b = 17°, the 
expected velocity range for fountain clouds is —30 to 15 
km s _1 (note that clouds at the positive velocities would 
be very difficult to detect due to confusion with Galactic 
HI), again very much in agreement with what we find for 
LVCs. The observed iV(HI) of LVCs is about ten times 
higher that what is expected for fountain clouds, although 
Houck & Brcgman (1990) note that A^HI) could be higher 
if clouds are concentrated at lower z heights. 



Another important evidence in favor of the fountain- 
style clouds is hinted by cloud deviation velocities. Collins 
et al. (2002) investigated a ballistic infall of fountain clouds 
and showed that a vertical velocity gradient (with scale 
height z) is expected for halo clouds; this is a consequence 
of a decrease in the gravitational acceleration vector with 
height from the disk. Hence, at any galactocentric ra- 
dius, the LSR velocity of clouds in the Halo, at z > 1 
kpc, will be different from the LSR velocity of clouds in 
the disk. However, the steepness of this vertical veloc- 
ity gradient is a function of the galactocentric radius at 
which clouds originated before being kicked upward from 
the disk. For example, for the input circular velocity of 
230 km s _1 and at a galactocentric radius of R g — 8 kpc, 
the expected AVi sr is <J 30 km s _1 when z changes from 
z = to z <~ 2 kpc (assuming an initial cloud kick velocity 
of 100 km s -1 ). Note that viscous interactions between 
Halo clouds, and/or more complex interactions between 
the Halo medium and clouds can decrease AVj sr . However, 
for very small galactocentric radii, AVj sr « km s _1 , for 
any z. The expected LSR velocity difference (AVi sr ) for 
clouds at higher z relative to clouds in the disk, is there- 
fore something we can directly compare with the cloud 
deviation velocity, Vdov 

The anti-center LVCs have (|Vdov|) = 30 km s _1 , in 
agreement with what is expected for fountain clouds at 
heights z =1-2 kpc. The expected circulation timcscale 
(time for clouds to return back to the disk) for the fountain 
clouds is about 50 Myr (Houck & Bregman 1990; Collins 
et al. 2002). As the cloud evaporation timescale (;> 100 
Myr) is significantly longer than the expected circulation 
timescale, clouds can reach the disk safely before being 
evaporated. In Section 7.2 we hinted at a potentially in- 
teresting difference in |Vd ev | for clouds in the inner and 
outer Galaxy — most halo clouds in the inner Galaxy ap- 
pear to have smaller |Vdev| than clouds in the outer Galaxy. 
This potential trend in |Vdcv| is an important tool for in- 
vestigating cloud origin — if halo clouds are a result of a 
Galactic fountain we would expect |Vd cv | to increase with 
galactocentric radius. Future work with larger cloud sam- 
ples will investigate statistical significance of this trend. 

In conclusion, LVCs have typical linear sizes and ra- 
dial velocities in good agreement with expectations from 
the Galactic fountain model. This would suggest that 
Galactic fountain flows, and the circulation of fountain 
clouds, are (at least) as common in the outer Galaxy as 
in the inner Galaxy. This is a surprising result when 
one takes into consideration that the Galactic star forma- 
tion rate (SFR) strongly depends on galactocentric radius, 
based on various observations as well as predictions from 
chemical and spectrophotometric models for the evolution 
of spiral galaxies (for references see Boissier & Prantzos 
(1999) and Boissier et al. (2003)). The star formation 
rate peaks around R g ~ 4 kpc and then steadily de- 
creases towards larger radii, at R g ~ 12 kpc for example, 
SFR— 0.3 x SFR Q . However, the decrease in SFR at large 
radii is alleviated partially by the decrease of the gravi- 
tational potential which enables vertical fountain flows to 
propagate easier. 

9.2. Connection between filaments and clouds 
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Most of LVCs appear related to large, diffuse HI fila- 
mentary features. Diffuse filaments and clouds both ap- 
pear only over a restricted range of LSR velocities, —20 to 
—30 km s _1 . This is suggestive of their possible physical 
association. 

An association between larger filaments and small, cold 
clouds is expected as a result of shell or chimney frag- 
mentation. Norman & Ikeuchi (1989) expect that clouds 
formed from chimney walls will have velocities \v\ ~10-50 
km s" 1 . McClure-Griffiths et al. (2006) found several cold 
HI clouds associated with caps of a large expanding shell 
and suggested that shell fragmentation is responsible for 
the cloud formation. Similarly, Koo et al. (2006) proposed 
that small, fast-moving clouds could be fragments of old 
supernova remnant shells. 

Visual appearance of filaments and clouds (Figure 4) 
strikingly resembles recent numerical simulations by Au- 
dit & Henncbclle (2005) of dynamically triggered conden- 
sation of the WNM into small CNM clouds. In these 
simulations, a collision between incoming turbulent WNM 
streams creates a thermally unstable region of higher den- 
sity and pressure but lower temperature, which starts to 
fragment into cold structures. This thermally unstable 
WNM has a filamentary morphology, and is seeded with 
numerous small (size of ~ 0.1 pc) and cold clouds (n ~ 50 
cm~ 3 , T ~ 80 K). The whole process of fragmentation is 
promoted and controlled by turbulence, and cloud mor- 
phology and properties, as well as the fraction of the con- 
densed cold gas, are governed by turbulent properties. 
Based on this apparent morphological similarity it may 
be reasonable to expect that external (dynamical) forc- 
ing (through, for example, collisions, or interactions with 
the dense halo gas) of some types of vertical WNM flows 
(e.g. chimneys or walls of expanding shells) can trigger 
condensation of warm HI into cold clouds. 

9.3. Inf ailing multi-phase gas? 

The idea that galaxies acquire their fresh star forma- 
tion fuel through the infall of clouds has attracted a lot of 
attention, especially recently with advances in numerical 
simulations of galaxy formation (Mailer & Bullock 2004). 
As the hot gas within the cooling radius (~ 150 kpc for the 
Milky Way size galaxy) is a subject to thermal instabilities, 
it fragments into smaller, warm (T ~ 10 4 K), pressure- 
supported clouds. The warm clouds fall into the Galaxy, 
after undergoing cloud-cloud collisions and ram pressure 
stripping, and are tidally broken into smaller clouds before 
impacting the Galaxy. The cloud formation and infall are 
expected to be balanced at the present time, and models 
predict abundant evidence for this ongoing process in the 
Galactic halo, in the form of numerous condensed clouds. 

Kaufmann et al. (2006) investigate the formation of a 
galactic disk, through N-body and SPH simulations, and 
show that cloud accretion is not spherical but happens 
along the angular momentum axis. This results in outer 
disks preferentially being populated with small, cold accre- 
tion remnants, while the central regions below and above 
the plane channeling hot gas. In this scenario, cloudy 
structure of the disk/halo interface region is expected in 
the inner and outer galactic disks but due to a very dif- 
ferent phenomenon — while outer parts are driven by the 
accretion process, the inner parts could be a result of a 



fountain process. 

However, what exactly happens to clouds during the fi- 
nal stages of their infall is still mainly unexplored, both 
theoretically and observationally. How do large and warm 
clouds transform into 'cloud forms' that can easily be in- 
tegrated and digested by the disk, to provide a smooth 
build up of fresh star formation fuel? Mailer & Bullock 
(2004) suggest that clouds get shredded into smaller but 
still warm (~ 10 4 K) condensations, with further cooling 
being stopped by the extragalactic ionizing background. 
Clouds in outer disks in simulations by Kaufmann et al. 
(2006) are also warm, 100-600 pc in size, however appar- 
ently they can cool and condense further to reach thermal 
equilibrium at a temperature < 10 4 K. 

CHVCs and our LVCs may be able to fit into this 
framework as representing later stages of the cloud in- 
fall process. CHVC186+19-114 shown in Figures 8-11, 
typical of several hundreds of compact HVCs, could rep- 
resent a cloud that has condensed in the halo and has 
begun to infall, but has not yet reached the disk. The 
low- velocity filamentary/cloudy structure, shown in Fig- 
ures 2-7, could represent the remaining structure of the 
accreted cloud as it is being almost integrated with the 
Galactic disk. This scenario seems plausible when one con- 
siders the similarities between these two types of objects 
(see Table 1): CHVC186+19-114's 'core' structure has ve- 
locity linewidths and HI column densities comparable to 
those found in LVCs, while typical HI linewidths and col- 
umn densities of the filamentary structures at about —25 
km s" 1 are very similar to those of CHVC186+19-114's 
envelope. If the LVCs trace the final stages of the infall 
process then they suggest cloud sizes and temperatures 
significantly smaller from what is currently found in nu- 
merical simulations. However, details of heating/cooling 
and self-shielding effects during the cloud infall still await 
to be addressed. Observationally, constraining metallici- 
ties of halo clouds would be the crucial test of the cloud 
infall hypothesis. 

10. SUMMARY AND FUTURE WORK 

The consortium for Galactic studies with ALFA is con- 
ducting an HI survey of the whole Arecibo sky (declination 
range from —1° to 38°) with high angular (3.5') and ve- 
locity resolution (0.2 km s _1 ). At the time this article is 
being written about 20% of the survey has already been 
completed. In this paper we focused only on one of the 
main science drivers for GALFA — the Galactic disk/halo 
interface region. With only about 200 square degrees im- 
aged in the Galactic anti-center region, primarily for the 
purpose of hardware and software testing, several interest- 
ing and new phenomena have emerged. 

We have found numerous small (most likely just a few pc 
in size) and cold (Tk < 400 K) HI clouds at low negative 
velocities distinctly separated from the HI disk emission 
('low- velocity clouds' or LVCs). Over a restricted range 
of velocities and exactly the same range as for LVCs, we 
have also found larger sheet-like filamentary structures. 
Most likely, the LVCs and filaments are physically related. 
While the distances to LVCs will have to be better con- 
strained in the future, preliminary estimates, based on 
the thermal pressure requirements for the CNM and the 
WNM coexistence, range from 0.2 to 3 kpc. This implies 
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that the clouds' height from the plane is in the range 60 
to 900 pc, suggesting that most LVCs are located in the 
transition region between the Galactic disk and halo, yet 
they have properties of typical CNM clouds. LVCs are 
colder and, most likely, smaller and less massive than Lock- 
man's clouds in the disk/halo interface region in the inner 
Galaxy. In addition, LVCs may have larger absolute devi- 
ation velocities than Lockman's clouds. Nevertheless, the 
existence of a large number of LVCs in the outer Galaxy 
(at a galactocentric radius > 8 kpc) demonstrates that the 
cloudy and frothy character of the interface region is more 
likely to be a general phenomenon, not restricted exclu- 
sively to the inner Galaxy. 

In the same dataset but at more extreme negative 
velocities we have studied CHVC186+19-114, a typical 
CHVC with a well-defined core/envelope structure. As 
we have imaged a large area around the CHVC, we 
discovered a companion HI cloud located 50' southwest 
of CHVC186+19-114, most likely physically related to 
CHVC186+19-114. Velocity profiles suggest that the com- 
panion cloud was stripped off the main cloud's envelope 
due to the interactions with the halo medium. This is 
the first time that a companion is detected around an iso- 
lated CHVC and may represent evidence for the breakup 
of CHVCs into smaller condensations. Based on its size, 
the companion cloud could be classified as a member of 
the population of UCHVCs, discovered only very recently. 

We have discussed briefly three possible scenarios for 
the formation and maintenance of the clumpy Galactic 
disk/halo interface region. The first scenario is provided 
by Galactic fountains, buoyant hot flows from the disk 
which cool and rain back in the form of cold clouds. Ob- 
servational properties (cloud size and radial velocity) of 
LVCs are in a good agreement with predictions of low- 
temperature fountains models. The second scenario is mo- 
tivated by numerical simulations by Audit & Hennebelle 
(2005) of dynamically triggered condensation of the WNM 
into small CNM clouds. Vertical WNM flows propagating 
upward from the Galactic plane, in the form of expanding 
shells or chimneys, could be triggered into formation of 
CNM clouds at larger z heights. This comparison is pri- 
marily based on morphological resemblance at this stage as 
detailed models do not exist. Finally, the third scenario in- 
vestigates whether the low- velocity clouds/filaments could 
represent the final stages of the infalling IGM. This has 



been encouraged by the recent models and simulations 
(Mailer & Bullock 2004; Kaufmann et al. 2006) which em- 
phasize the multi-phase and multi-scale character of the 
galaxy formation process. While simulations show encour- 
aging qualitative resemblance at this stage the detailed 
physical processes during the final stages of cloud infall 
are largely unexplored. We hope that our current, and es- 
pecially future high- resolution datasets, will stimulate and 
guide further theoretical and numerical work in this area. 

Future GALFA datasets will clearly find more disk/halo 
interface clouds and new HVCs, especially in the regime 
of the smallest HVCs which is largely unexplored. With 
the larger data sets we will be able to confirm whether the 
cloudy structure exists at other Galactic longitudes, and 
also the statistical significance of the potential differences 
between interface clouds in the inner and outer Galaxy. 
Also, future GALFA datasets will be important for estab- 
lishing spatial distribution and properties of the smallest 
HVCs. 
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